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Abstract 

The NMSSM contains a Higgs singlet in addition to the two Higgs doublets typical of 
the MSSM, thus resulting in a total of seven physical Higgs mass states. Therefore, 
the phenomenology of the NMSSM Higgs sector can vary considerably from that of the 
MSSM and there are good prospects of finding in regions of the NMSSM parameter 
space Higgs signals that cannot be reproduced in the MSSM. We examined here the 
two-photon decay mode of a Higgs boson and found that up to three neutral Higgs 
states, heavy and/or light, could be simultaneously observable at the LHC, a possibility 
precluded to the MSSM. There are also some possibilities that only the lightest NMSSM 
Higgs boson be detectable via this mode, with a mass beyond the upper limit of the 
corresponding MSSM state, thus also allowing to distinguish between the two scenarios. 
However, in most of the NMSSM parameter space the configurations of the non-minimal 
model are not very different from those arising in the minimal case. 



1 Introduction 

The Minimal Supersymmetric Standard Model (MSSM) lj is affected by the so-called '/i- 
problem'. The MSSM Superpotential can be written as j2] 

^mssm = QH u h u U c + H d Qh d D c + H d Lh e E c + fiH u H d . (1) 

(Hereafter, hatted variables describe Superflelds while un-hatted ones stand for the corre- 
sponding scalar Superfield components). The last term in the above equation contains a 
dimensionful parameter, /x. Upon Electro- Weak Symmetry Breaking (EWSB), it provides a 
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contribution to the masses of both Higgs bosons and Higgsino fermions. Furthermore, the 
associated soft Supersymmetry (SUSY) breaking term BfiH u Hd mixes the two Higgs dou- 
blets. Now, the presence of /i in the Superpotential before EWSB indicates that its natural 
value would be either or the Planck mass Mp. On the one hand, \i = would mean 
no mixing is actually generated between Higgs doublets at any scale and the minimum of 
the Higgs potential occurs for < Hd >= 0, so that one would have in turn massless down- 
type fermions and leptons after SU(2) symmetry breaking. On the other hand, /i ~ Mp 
would reintroduce a 'fine-tuning problem' in the MSSM since the Higgs scalars would ac- 
quire a huge contribution ~ jj 2 to their squared masses (thus spoiling the effects of SUSY, 
which effectively removes otherwise quadratically divergent contributions to the Higgs mass 
from SM particles). Therefore, the values of this (theoretically arbitrary) parameter fi are 
phenomenologically constrained to be close to Msusy or Mw 

The most elegant solution to the /x-problem is to introduce a new singlet Higgs field S 
and replace B\iH u H d by an interaction ~ S(H u H d ). When the extra scalar field S acquires a 
Vacuum Expectation Value (VEV), an effective /i term, naturally of the EW scale, is gener- 
ated automatically. This idea has been implemented in the Next-to-Minimal Supersymmetric 
Standard Model (NMSSM) described by the Superpotential 

Wnmssm = QH u h u U c + H d Qh d D c + H d Lh e E c + XS(H u H d ) + ^kS\ (2) 

where S is an extra Higgs iso-singlet Superfield, A and k are dimensionless couplings and 
the last (Z 3 invariant) term is required to explicitly break the dangerous Peccei-Quinn (PQ) 
U(l) symmetry [Sj 1 . (See Ref. [Zj for NMSSM Higgs sector phenomenology with an exact or 
slightly broken PQ symmetry.) Furthermore, a symmetry can be imposed to avoid the 
so-called 'domain- walls problem' of the NMSSM [H], through harmless tadpoles breaking the 
global Z 3 symmetry |3J. (Alternative formulations - known as the Minimal Non-minimal 
Supersymmetric Standard Model (MNSSM) and new Minimally-extended Supersymmetric 
Standard Model (nMSSM) - vetoing the presence of the ~ kS 3 term and allowing instead 
suitable linear tadpole terms through the enforcement of discrete i?-symmetries - also exist 
|10j.) Another positive feature of all these non-minimal SUSY models is that they predict the 
existence of a (quasi-) stable singlet-type neutralino (the singlino) that could be responsible 
for the Dark Matter (DM) of the universe [IT] . Finally, notice that, in these extended 
SUSY models, the singlet Superfield S has no SM gauge group charge (so that MSSM 
gauge coupling unification is preserved) and that one can comfortably explain the baryon 
asymmetry of the Universe by means of a strong first order EW phase transition ^2] (unlike 
the MSSM, which would require a light top squark and Higgs boson barely compatible with 
current experimental bounds fTHj ) . 

Clearly, in eq. (0), upon EWSB a VEV will be generated for the real scalar component of 
5* (the singlet Higgs field), < S >, alongside those of the two doublets < H u > and < Hd > 
(related by the parameter tan/3 =< H u > / < Hd >). In the absence of fine-tuning, one 
should expect these three VEVs to be of the order of M SU sy or M w , so that now one has an 

1 One could also gauge the U(1)pq group, so that the Z 3 symmetry is embedded in the local gauge 
symmetry 0. 
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'effective /^-parameter' 

yUcff = A < S >, 
of the required size, thus effectively solving the /z-problem. 



(3) 



2 Higgs Phenomenology in the NMSSM 

In the NMSSM, the soft SUSY-breaking Higgs sector is described by the Lagrangian contri- 
bution 

Vnmssm = m 2 H jH u \ 2 + m 2 Hd \H d \ 2 + m 2 s \S\ 2 + (\A x SH u H d + ^A K S 3 + h.c.J , (4) 

where A\ and A K are dimensionful parameters of order Msusy- 

As a result of the introduction of an extra complex singlet scalar field, which only couples 
to the two MSSM-type Higgs doublets, the Higgs sector of the NMSSM comprises of a total of 
seven mass eigenstates: a charged pair h ± , three CP-even Higgses ^1,2,3 ( m hi < m h 2 < m h z ) 
and two CP-odd Higgses a 12 (m ai < m a2 ). Consequently, Higgs phenomenology in the 
NMSSM may be plausibly different from that of the MSSM and extremely rich of new 
signals. 

For a start, the mass expressions for the CP-even Higgses in the NMSSM can be translated 
into a modified upper bound on the lightest Higgs mass, m hl , as [H] 

m 2 hl < min{m|, -k < S > (Ak < S > +V2A K )}. (5) 

As the higher order corrections are similar to those in the MSSM, the upper bound on the 
lightest Higgs boson is different in the NMSSM, reaching 135-140 GeV, for maximal stop 
mixing and tan/? = 2 ^3 IE] (a configuration indeed excluded in the MSSM by LEP data). 
More in general, the 'little fine tuning problem', resulting in LEP failing to detect a light CP- 
even Higgs boson, predicted over most of the MSSM parameter space, is much attenuated 
in the NMSSM, because the mixing among more numerous CP-even or CP-odd Higgs fields 
enables light mass states being produced at LEP yet they can remain undetected because of 
their reduced couplings to Z bosons [THj . 

As for future machines, chiefly the CERN Large Hadron Collider (LHC), quite some 
work has been dedicated to probing the NMSSM Higgs sector over recent years. Primarily, 
there have been attempts to extend the so-called 'No-lose theorem' of the MSSM - stating 
that at least one MSSM Higgs boson should be observed via the usual SM-like production 
and decay channels at the LHC throughout the entire MSSM parameter space PG] - to the 
case of the NMSSM ^H]- From this perspective, it was realised that at least one NMSSM 
Higgs boson should remain observable at the LHC over the NMSSM parameter space that 
does not allow any Higgs-to-Higgs decay. However, when the only light non-singlet (and, 
therefore, potentially visible) CP-even Higgs boson, hi or h 2 , decays mainly to two very light 
CP-odd Higgs bosons, h — > a^i, one may not have a Higgs signal of statistical significance 
at the LHC ^Hj- In fact, further violations to the theorem may well occur if one enables 
Higgs-to-SUSY decays (e.g., into neutralino pairs, yielding invisible Higgs signals). 
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While the jury is still out on whether a 'No-lose theorem' can be proved for the NMSSM, 
we are here concerned with an orthogonal approach. We asked ourselves if a, so to say, 'More- 
to-gain theorem' can be formulated in the NMSSM. That is, whether there exist regions of 
the NMSSM parameter space where more Higgs states of the NMSSM are visible at the LHC 
than those available within the MSSM. In our attempt to overview all such possibilities, we 
start by considering here the case of the di-photon decay channel of a neutral Higgs boson. 
This mode can be successfully probed in the MSSM, but limitedly to the case of one Higgs 
boson only, which is CP-even and rather light. We will show that in the NMSSM one can 
instead potentially have up to three di-photon signals of Higgs bosons, involving not only 
CP-even but also CP-odd states, the latter with masses up to 600 GeV or so. In fact, even 
when only one di-photon signal can be extracted in the NMSSM, this may well be other than 
the h\ state. When only the latter is visible, finally, it can happen that its mass is larger 
than the maximum value achievable within the MSSM. In all such cases then, the existence 
of a non-minimal SUSY Higgs sector would be manifest. 

3 Parameter Space Scan 

The choice of parameter space largely depends on the version of NMSSM under considera- 
tion, or more specifically, on the implication of unification of parameters at some very high 
scale. This, in turn, leans on the technique adopted for the breaking of SUSY, since it 
still remains undetected in nature. The advantage of assuming unification of masses and 
couplings in a SUSY model is twofold. Firstly, it enormously reduces the labour from a phe- 
nomenological perspective by minimising the number of parameters required to extrapolate 
physical information from the model. Secondly, and theoretically more crucially, it caters to 
the fundamental objective of constructing a Grand Unification Theory (GUT), which was, 
to a great extent, responsible for devising SUSY in the first place. 

A particular case of a low energy NMSSM with gravity mediated SUSY breaking was 
studied in [20] • As usual, in such a model, SUSY is assumed to be broken in some hidden 
sector and then mediated to the physical sector through gravitational interactions. Such a 
model implies unification of the couplings and soft masses at the GUT scale, which are then 
run down, using renormalisation group equations, to the weak scale or to some other scale at 
which the theory is being tested. This results in only a handful of parameters to deal with. 

For a more general study of the NMSSM Higgs sector with a wider range of parameters, 
we used here the publicly available fortran code NMHDECAY (version 1.1) developed in 
Ref. [21 J. This program computes the masses, couplings and decay widths of all the Higgs 
bosons of the NMSSM in terms of its parameters at the EW scale. For our purpose, instead 
of postulating unification, and without taking into account the SUSY breaking mechanism, 
we fixed the soft SUSY breaking terms to a very high value, so that they have little or no 
contribution to the outputs of the parameter scans. Consequently, we are left with six free 
parameters. 

Our parameter space includes the Yukawa couplings A and k, the soft trilinear terms A\ 
and A K , plus tan/3 and /x e g = X(S). The computation of the spectrum includes leading two 
loop terms, EW corrections and propagator corrections. The decay widths, however, do not 
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include three body decays. The NMHDECAY program also takes into account theoretical 
as well as experimental constraints from negative Higgs searches at LEP, along with the 
unconventional channels relevant for the NMSSM. 

We have used the NMHDECAY code to scan over the NMSSM parameter space defined 
through the aforementioned six parameters taken in the following intervals: 

A : 0.0001 - 0.75, « : -0.65 - +0.65, tan/5 : 1.6 - 54, 
//, A x , A K : -1000 - +1000 GeV. 

Remaining soft terms which are fixed in the scan include: 

• m Q3 = m Ua = m D:j = m La = m Ea = 2 TeV, 

• A U3 = A D3 = A E3 = 1.5 TeV, 

• mg = mu = rriD = mi = mE = 2 TeV, 

• Mi = M 2 = M 3 = 3 TeV. 

In line with the assumptions made in [22], the allowed decay modes for neutral NMSSM 
Higgs bosons are 2 : 

h,a^gg, h,a—*[i + fi~, h,a-^T + r~, h,a-^bb, h,a—*ti, 
h,a —>■ ss, h,a — > cc, h — > W + W~ , h — > ZZ, 
h, a — > 77, h, a — > Zj, h,a — > Higgses, h,a sparticles. 

(Notice that for the pseudoscalar Higgses, the decay into vector boson pairs is not allowed 
due to CP-conservation.) Here, 'Higgses' refers to any final state involving all possible 
combination of two Higgs bosons (neutral and/or charged) or of one Higgs boson and a 
gauge vector. 

The aforementioned range of parameters is borrowed from Ref. [22], where the authors 
singled out the most difficult scenarios for NMSSM Higgs discovery at the LHC (using the 
mentioned code). In this paper they also concluded, as a follow up of the NMSSM 'No- 
Lose theorem', that the region where no Higgs is observable, due to its decays into the 
unconventional channel h — > a^x, comprises of only about 1% of the entire parameter space 
of the NMSSM. In short, our idea is to explore the same range of parameter space and look 
for the 'best case scenarios' where the discovery of one or more Higgses at the LHC is possible 
above and beyond what is predicted in the MSSM. 

We have performed our scan over several millions of randomly selected points in the 
specified parameter space. The output, as stated earlier, contains masses, Branching Ratios 
(BRs) and couplings of the NMSSM Higgses, for all the points which are not forbidden by the 
various constraints. The points which violate the constraints are eliminated. The surviving 
data points are then used to determine the cross-sections for NMSSM Higgs hadro-production 
by using an adapted version of the codes described in [23] ■ As the SUSY mass scales have 
been arbitrarily set well above the EW one (see above), the NMSSM Higgs production modes 
exploitable in simulations at the LHC are those involving couplings to heavy ordinary matter 

2 Here, we use the label h(a) to signify any of the neutral CP-even(odd) Higgs bosons of the NMSSM. 
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only, i.e., (hereafter, V = W ± , Z and Q = b,t) for neutral Higgs production (where the last 
two channels are only allowed for CP-even Higgs production): 

gg — > Higgs (gluon — fusion, via heavy — quark loops), 

99 ~ > QQ Higgs (heavy — quark associated production), 

qq — > qqV*V* — > qq Higgs (vector — boson — fusion), 

qq — > V Higgs (Higgs — strahlung). 

(These are the so-called 'direct' Higgs production modes.) Here, 'Higgs' refers to any possible 
neutral Higgs boson. 

Production and decay rates for NMSSM neutral Higgs bosons have then been multiplied 
together to yield inclusive event rates, assuming a LHC luminosity of 100 fb _1 throughout. 

4 Spectrum Configuration and Inclusive Event Rates 

As an initial step towards the analysis of the data, we computed the total cross-section 
times BR into 77 pairs against each of the six parameters of the NMSSM, for each Higgs 
boson. We have assumed all production modes described in the above Section and started 
by computing total (i.e., fully inclusive) rates 3 . We are here focusing on the 77 decay mode 
since it is the most promising channel for the discovery of a (neutral) Higgs boson at the 
LHC in the moderate Higgs mass range (say, below 130 GeV). In fact, since the tail of 
the 77 background falls rapidly with increasing invariant mass of the di-photon pair, signal 
peaks for heavier Higgses could also be visible in addition to (or instead of) the lightest one, 
although the cross-section for these processes is relatively very low |2H EE] • As the starting 
point of our numerical study, based on the ATLAS analysis of Ref . [2H] , we argue that cross- 
section times BR rates of 10 fb or so are potentially interesting from a phenomenological 
point of view, in the sense that they may yield visible signal events, the more so the heavier 
the decaying Higgs state (also because the photon detection efficiency grows with the Higgs 
mass [2E])- (See Ref. |2Zj for a preliminary account in this respect.) 

Tab. ^ recaps the potential observability of one or more NMSSM Higgs states in the 
di-photon mode at the LHC, under the above assumptions. It is obvious from the table that 
one light CP-even Higgs should be observable almost throughout the NMSSM parameter 
space (in line with the findings of Ref. [22] )• However, there is also a fair number of points 
where two Higgses may be visible simultaneously (hi and hi or - more rarely - h\ and 
ai), while production and decay of the three lightest Higgses (h%, h 2 and ai) at the same 
time, although possible, occurs for only a very small number of points in the parameter 
space. Furthermore, the percentage of points for which only the second lightest Higgs state 
is visible is also non- negligible. These last two conditions are clearly specific to the NMSSM, 
as they are never realised in the MSSM. Furthermore, while the lone detection of the lightest 
CP-even NMSSM Higgs boson may mimic a similar signal from the corresponding state in 

3 After verifying that the bulk of the signal rates is due to gluon-gluon fusion (even at large Higgs masses), 
we have eventually decided - for simplicity - to limit ourselves to emulate only this channel. Hence, all the 
results found below suffer from a slight under- estimate of the signal rates. 
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the MSSM, the reconstructed mass may well be beyond the upper mass limit in the MSSM, 
this possibility also pointing towards the evidence of a NMSSM Higgs sector. Finally, none 
of the two heaviest NMSSM neutral Higgs states (h 3 and a 2 ) will be visible in the di-photon 
channel at the LHC (given their large masses). 

We have then plotted the cross-sections times BR for the three potentially observable 
Higgses against the various NMSSM parameters. These plots, shown in Figs. 1-6, reveal that 
the distribution of possibly visible points (i.e., of those yielding cross-section time BR rates 
in excess of 10 fb) is quite homogeneous over the NMSSM parameter space and not located in 
some specific parameter areas (i.e., in a sense, not 'fine-tuned'). The distribution of the same 
points as a function of the corresponding Higgs masses can be found in Fig. 7. Of particular 
relevance is the distribution of points in which only the NMSSM hi state is visible, when 
its mass is beyond the upper mass limit for the corresponding CP-even MSSM Higgs state, 
which is shown in Fig. 8 4 . This plot reveals that about 93% of the NMSSM h\ masses visible 
alone are expected to be within 2-3 GeV beyond the MSSM bound, hence the two models 
would be indistinguishable 5 . Nonetheless, there is a fraction of a percent of such points 
with mh 1 values even beyond 125 GeV or so (the higher the mass the smaller the density, 
though), which should indeed allow one to distinguish between the two models. Moreover, by 
studying the cross-section times BR of the Higgses when two of them are observable against 
their respective mass differences, Figs. 9-11, one sees that the former are larger than the 
typical mass resolution in the di-photon channel, so that the two decaying objects should 
indeed appear in the data as separate resonances. (We have also verified, though not shown 
here, that their decay widths are small compared to the detector resolution in M 77 .) 

Next, we have proceeded to a dedicated parton level analysis of signal and background 
processes, the latter involving both tree-level qq — > 77 and one-loop gg — > 77 contributions. 
We have adopted standard cuts on the two photons [2EJ: j? T > 25 GeV and |?7 7 | < 2.4 on 
transverse energy and pseudorapidity, respectively. As illustrative examples of a possible 
NMSSM Higgs phenomenology appearing at the LHC in the di-photon channel, we have 
picked up the following three configurations: 

1. A = 0.6554, k = 0.2672, /i = -426.48 GeV, tan/3 = 2.68, A x = -963.30 GeV, 
A K = 30.48 GeV; 

2. A = 0.6445, k = 0.2714, // = -167.82 GeV, tan/3 = 2.62, A x = -391.16 GeV, 
A K = 50.02 GeV; 

3. A = 0.4865, k = 0.3516, /1 = 355.63 GeV, tan/3 = 2.35, A x = 519.72 GeV, A K = 
-445.71 GeV. 

"Notice that the value obtained for m^ ax from NMHDECAY version 1.1, of - 130 GeV, based on the 
leading two-loop approximations described in is a few GeV lower than the value declared in Sect. 1. 
Besides, for consistency, we use the value of 120 GeV (obtained at the same level of accuracy) as upper 
mass limit on the lightest CP-even Higgs boson of the MSSM. (Notice that a slightly modified m™ ax value is 
obtained for the NMSSM from NMHDECAY version 2.1 [2E], because of the improved mass approximations 
with respect to the earlier version of the program adopted here.) Eventually, when the LHC is on line, the 
exercise that we are proposing can be performed with the then state-of-the-art calculations. 

5 Other than an experimental di-photon mass resolution of 2 GeV or so one should also bear in mind 
here that the mass bounds in both models come with a theoretical error of comparable size. 



7 



The first is representative of the case in which only the NMSSM hi boson is visible, but 
with mass larger than the MSSM upper limit on the corresponding Higgs state. The second 
and third refer instead to the case when also the h 2 or a\ state are visible, respectively. 
The final results are found in Fig. 12. The corresponding mass resonances are clearly visible 
above the continuum di-photon background and discoverable beyond the 5a level. Indeed, 
similar situations can be found for each of the combinations listed in Tab. Q and most of 
these correspond to phenomenological scenarios which are distinctive of the NMSSM and 
that cannot be reproduced in the MSSM. 

5 Conclusion 

In summary, we have shown that there exists the possibility of establishing a 'More-to-gain 
theorem' within the NMSSM, as compared to what is expected in the MSSM, in terms of 
novel Higgs signals appearing in the di-photon discovery channel which can be ascribed to the 
former but not to the latter. We have showed this to be the case for a few selected NMSSM 
parameter points, by performing a proper signal-to-background analysis at the partonic level. 
However, a similar numerical study can easily be extended to encompass sizable regions of the 
NMSSM parameter space. While the bulk of the latter is in a configuration degenerate with 
the MSSM case (as far as di-photon signals at the LHC are concerned), non- negligible areas 
exist where further phenomenological studies have the potential to unveil a non-minimal 
nature of the underlying SUSY model. 

To this end, NMSSM benchmark scenarios, amenable to experimental investigation in 
the context of full Monte Carlo (MC) analyses, based on event generation and detector 
simulation, are currently being devised Also, un upgraded version of HERWIG [SO] , 

suitable for MC event generation in the context of non-minimal SUSY models, is currently 
being prepared [8T] . 
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Higgs Flavor 


Points Visible 


Percentage 




Total: 


1345884 


99.7468 




Alone: 


1345199 


99.6961 


h 


With h 2 : 


528 


0.0391 




With ay. 


152 


0.0113 




With h 2 and ay. 


5 


0.0004 




Total: 


1253 


0.0929 




Alone: 


717 


0.0531 


h-2 


With hy. 


528 


0.0391 




With ay. 


3 


0.0002 




With ai and ay. 


5 


0.0004 


hs 


Total: 










Total: 


165 


0.0122 




Alone: 


5 


0.0004 


ai 


With hy 


152 


0.0113 




With h 2 : 


3 


0.0002 




With hi and /i 2 : 


5 


0.0004 


a 2 


Total: 









Table 1: Higgs events potentially visible at the LHC through the 77 decay mode (i.e., those 
yielding cross-section time BR rates of order 10 fb or upwards). Percentage refers to the 
portion of NMSSM parameter space involved for each discovery scenario. 
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Figure 1: Cross-section times BR of /ii (red), /i 2 (green) and a x (blue), when potentially 
visible individually and when two of these are potentially visible simultaneously, plotted 
against the parameter A. 
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Figure 2: Cross-section times BR of hi (red), h% (green) and a\ (blue), when potentially 
visible individually and when two of these are potentially visible simultaneously, plotted 
against the parameter k. 



13 




-ICO -2SC 




—SCO — ICC -3CO -2CO -ICO ICO 2C0 3C0 4CO 




w [G.V] 



Figure 3: Cross-section times BR of h\ (red), /i 2 (green) and a\ (blue), when potentially 
visible individually and when two of these are potentially visible simultaneously, plotted 
against the parameter \i. 
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Figure 4: Cross-section times BR of hi (red), /i 2 (green) and a x (blue), when potentially 
visible individually and when two of these are potentially visible simultaneously, plotted 
against the parameter tan/3. 
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Figure 5: Cross-section times BR of hi (red), h 2 (green) and ai (blue), when potentially 
visible individually and when two of these are potentially visible simultaneously, plotted 
against the parameter A\. 
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Figure 6: Cross-section times BR of hi (red), h% (green) and a\ (blue), when potentially 
visible individually and when two of these are potentially visible simultaneously, plotted 
against the parameter A K . 
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Figure 7: Cross-section times BR of hi (red), fo 2 (green) and ai (blue), when potentially 
visible individually and when two of these are potentially visible simultaneously, plotted 
against their respective masses. 
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Figure 8: The distribution of points where one has only one potentially visible NMSSM hi 
state with mass beyond the MSSM upper mass limit on the corresponding Higgs state. The 
scale on the right represents a measure of density of the points. 
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Figure 9: Cross-section times BR of h\ (red) and h 2 (green) plotted against their mass 
differences when the two are potentially visible simultaneously. 
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Figure 10: Cross-section times BR of (red) and a x (blue) plotted against thier mass 
differences when the two are potentially visible simultaneously. 
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Figure 11: Cross-section times BR of h 2 (green) and a\ (blue) plotted against their mass 
differences when the two are potentially visible simultaneously. 
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Figure 12: The differential distribution in invariant mass of the di-photon pair after the 
cuts in pj, and rp mentioned in the text, for 100 fb _1 of luminosity, in the case of the 
background (solid) and the sum of signal and background (dashed), for the example points 
1.-3. described in the text (from top to bottom, in correspondence). 
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